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U-spin multiplet approach is applied to the full set of charmless hadronic — > M^M^ decays for 
the purpose of precise extraction of the unitarity angle 7. Each of the four data sets, P°P^ , P^V*, 
V^P^ and V''^V^ , with P = pseudoscalar and V = vector, can be used to yield a precise value of 7. 
The crucial advantage of this method over the common SU(3) symmetry based quark-diagrammatic 
approach is that no assumptions regarding relative sizes of topological decay amplitudes need to 
be made. As a result, this method avoids an uncontrollable theoretical uncertainty that is related 
to the neglect of some topological diagrams (e.g., exchange and annihilation graphs) in the SU(3) 
approach. Application of the U-spin approach to the current data yields: 7 = (54lJ^^) . We find 

that improved measurements of (f)n^ and K ^ branching ratios would lead to appreciably better 
extraction of 7. In this method, which is completely data driven, in a few years we should be able 
to obtain a model independent determination of 7 with an accuracy of 0(few degrees). 



Precise determinations of the angles of the unitarity 
triangle (UT) remains an important but difScult goal in 
Particle Physics. Though methods for direct determina- 
tions of all the angles are now known, we are still quite 
far away from having large enough sample of B's that are 
needed The main challenge in extracting the angles 
from the data is of course that weak decays take place 
in the presence of strong interactions (i.e. QCD) which 
in this energy regime has important, non-perturbative 
effects. Fortunately, QCD respects flavor symmetries. 
Use of these symmetries presents an important avenue 
to extract results, though often at the expense of some 
accuracy. In the context of the angle 7 of the UT, in 
fact SU(3) flavor symmetry has already been put to very 
good use 2, 3]. In this paper we propose to use U-spin 
for determining 7 from charmless decays. 

There are substantial differences between U-spin multi- 
plet approach and other phenomenological methods, such 
as SU(3) based approach, of understanding the current 
B decay data. 

• U-spin multiplet method has the significant advan- 
tage that, unlike SU(3) fits to charmless B decays, 
quark diagrammatic topological approach is not in- 
voked at all. Thus, we do not need to make any 
assumptions about the relative sizes of various con- 
tributing topological diagrams and so no amplitude 
need be neglected 0- 

• U-spin is a flavor symmetry similar to isospin. 
Since it is a subgroup of SU(3) flavor symmetry, 
we expect it to be at least as accurate as SU(3) 
and possibly better. 

• It is important to emphasize that noticeable flavor 
symmetry breaking effects in decay amplitudes do 
not necessarily lead to large uncertainties in 7 ex- 
traction. For instance, SU(3) breaking effects of 
about 20% that are related to the ratio of decay 
constants Jk and /^r, only lead to a small (2°, or 
3%) theoretical uncertainty in determination of 7 



from SU(3) fits |5|. The effects of the rj — i]' mix- 
ing on the theoretical error in 7 were also found 
to be smafl (< 1°) 0. Since in the U-spin ap- 
proach graphical topologies are not used, estimate 
of U-spin breaking effects on 7 extraction may be 
amenable to calculational frameworks such as QCD 
factorization, pQCD, soft coUinear effective theory 
(SCET), or QCD sum rules 0]. 

• The fact that the U-spin approach does not make 
use of graphical topologies, of course means that 
electroweak penguins are automatically fully con- 
tained within this approach. This is in sharp con- 
trast to the case of a extraction with the use of 
isospin. Therein electroweak penguins complicate 
the extraction of the CKM phase a and present 
a serious limitation. In this respect the U-spin 
approach has an advantage over even the isospin 
method. 

• Needless to say, the standard B DK methods 
of direct 7 extraction are theoretically the cleanest 
(error of 0(.1%) IJ) and should ultimately pro- 
vide the most accurate determination of 7. But 
this accuracy will only be attained after very large 
data samples become available, perhaps many years 
down the road. The U-spin approach that we are 
using here, on the other hand, can provide a fairly 
accurate value of 7 (error of 0(few percent)) with 
modest increase of luminosities. 

• Furthermore, since the U-spin approach automat- 
ically includes all penguin contributions, whereas 
the B — > DK method do not, a comparison of the 
values of 7 from the two methods provides a crucial 
test for new physics. 

We will show that there are four separate sets of two- 
body decays of charged B's each of which can give a value 
of 7. Existing data already allows determination of 7 
with an accuracy in the same ball park as other methods 
being used. We identify modes (^tt^ and K* K^) whose 
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improved experimental measurements should appreciably 
improve the accuracy on 7 with this method. In the era 
of the current B-factories, with the planned luminosities 
of a few ab~^, the method should allow us to determine 7 
with an accuracy of a few degrees. Furthermore, as better 
experimental information, at these luminosities, becomes 
available for all the relevant data sets, this method should 
give an understanding of its inherent systematic error. 

Let us very briefly recapitulate some elementary as- 
pects of U-spin H, 0, 0, ■ Recall that the U-spin sub- 
group of SU (3) is similar to the I-spin (isospin) subgroup 
except that the quark doublets with U — 1/2, U3 = ±1/2 
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i?+ is a U-spin singlet, while charged charmless mesons 
7r+(p+), K+{K*+) belong to U-spin doublets. Neutral 
mesons may get contributions from a U-spin triplet and 
two U-spin singlets. Strange neutral mesons K'^{K*'~') 
and their antiparticles are pure U-spin triplets. The mul- 
tiplet decompositions of other neutral mesons can be de- 
termined to be 



tt" = -i|l 0) + ^|0 0)8 , 

rj = ^|1 0) + ^|0 0)8 -i|0 0)1 , 
= + 0)8 + ^10 0)1 , 
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where two U-spin singlets are defined as 

|0 0)8 = -^{ss + dd ~ 2uu) , 
73I " ■ 



|0 0)1 = -^\uu + dd + ss) 



(3) 



One may decompose the \AS\ = 1 and AS* = effec- 
tive Hamiltonians into members of the same two U-spin 
doublets multiplying given CKM factors. For practical 
purposes, using CKM unitarity, it is convenient to write 
the effective Hamiltonian so that it involves only the u 
and c quarks: 



In isospin analysis of B decays ^3 the effective Hamil- 
tonian transforms as either A/ = i or A/ — |. While 
electroweak penguins, violate isospin due to the charge 
difference between u and d quarks, they do not violate 
U-spin. There are only three topological diagrams that 
may contribute to charged B decays: tree, penguin (QCD 
and electroweak), and annihilation. The effective Hamil- 
tonian of any of these decay types transforms as a U-spin 
doublet, AU — ^. This makes U-spin a particularly con- 
venient approach that allows the complete description 
of charged B decays without making additional assump- 
tions on the size of individual topological diagrams and 
without neglecting any of them, including annihilation. 
While the SU(3) based approach 0,0] does not inherently 
require ignoring annihilation, exchange and penguin an- 
nihilation contributions, in practice one has to do that 
to limit the number of parameters and keep SU(3) fits 
stable. This advantage of the U-spin approach makes it 
particularly appealing; in the long run, it should signif- 
icantly reduce theoretical uncertainties associated with 
this method. 

Since the initial B^ meson is a U-spin singlet and 
the effective Hamiltonian always transforms like a U- 
spin doublet, the final M°M+ states must be U-spin 
doublets. They can be formed in three different ways. 
While the charmless charged meson M+ can only be- 
long to a doublet, the neutral meson M*^ can be a linear 
combination of three different multiplets. As a result, 
any AS" ~ 0, B^ M°M+ decay amplitude can be ex- 
pressed in terms of three amplitudes: Af, Aq, Bq. They 
correspond to the U-spin triplet, U-spin singlet |0 0)3, 
and SU(3) singlet |0 0)i contributions into the decay 
amplitude. Each of these three amplitudes consists of 
a "u-like" and a "c-like" part (Eq. EJ. Similarly, any 
I AS" I = 1 decay amplitude can be written in terms of 
three other amplitudes: Af, Aq, Bq. The assumption of 
U-spin symmetry implies that the difference between Af 
and Af comes only through the difference in the CKM 
matrix elements. Thus, the complete amplitudes for U- 
spin final states are given by 



A5 = : nl^'^ = V:,VudO-^ + V:,VM , (4) 
|A5| = 1 : nl^' = K;K.0." + V:,V,s01 . (5) 

The assumption of U-spin symmetry implies that U-spin 
doublet operators and are identical, as well as the 
O'ji and 01 operators. The subscripts d and s may be 
omitted. Hadronic matrix elements of these two opera- 
tors, and , will be denoted and A^ and will be 
referred to as "u-like" and "c-like" amplitudes , where 
the latter includes electroweak penguin contributions. 
Note that these amplitudes multiply different CKM fac- 
tors in |A5| = 1 and AS* = processes. 



A^ = : = V:^^aAl, + T^^KdAg,, 

\AS\ = 1 : Al, = V:,V,.sAl, + V;,V,,Al, 
B^ = V:,VusBS + V:,V,sBS . 



Then we find that physical decay amplitudes for V 
and y^y"*" modes may be decomposed into U-spin am- 
plitudes [Tl| . 
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A{K*°K+) 
A(cj7r+), 



A(K*°K* 
A{p"p+) 
A{ujp+) 
A{<t>p+) 



-) = 



A{K*° 



A{ujK*+^ 



-2V2Ai , 

-A^ -Ai + V2B^ 
V2Ad, + V2Ai + 
-2V2AI , 
3A^o + Al , 
-A^ + Al + V2B^ 



(7) 



Ai(l,K*+) = V2Af, - V2AI + B^ 



where Ai, Aq and Bq correspond to final states with vec- 
tor mesons in the U-spin triplet, in the octet U-spin 
singlet and in the SU(3) singlet, respectively. Naturally, 
the formulae for related V^P~^ and F'^F"'" decay modes 
are the same, as seen in the above relations. However, the 
actual values for each of the U-spin amplitudes are con- 
stant only within each of the two subsets. They accept 
different values in V'^P'^ and V'^V'^ subsets. 

Thus, eight V'^P^ decays are described by 12 parame- 
ters: six U-spin amplitudes |^o,il' l^o,il' l^o I ^^'^ l^ol' 
five relative strong phases between them and the weak 
phase 7. The same statement is separately valid for eight 
YOy+ ixLodes, too. 

In the same way one can decompose physical ampli- 
tudes for P'^P'^ and P'^V^ decay modes into U-spin am- 
plitudes. Then we find H^: 



A{K K+), 
A{Tr°n+), 
A{iTK+), 
A{v'n+), 

A{K°TT+), 

A(7rOi^+), 

A{vK+), 

Aiv'K+), 



A{K°K*- 
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Just as the two subsets of M'^AI'^ that were considered 
before, P^P^ and P^V^ are also separately described by 
similar 12 parameters. 

Charmless hadronic decays of the B'^ meson to the 
two-meson final states that contain vector V or pseu- 
doscalar P mesons comprise four subsets: P^P'^, V^V'^, 
F^P"*", and P°F+. Each of the subsets comprises eight 
decays, with all possible combinations of two charged 
mesons (e.g., 7r+ and K'^ in the pseudoscalar octet) and 
four neutral ones {K*°, p, lo, and (j) in the vector octet). 
Thus there are altogether 16 relevant decays of B^ of 
each of the four types. Each of the subsets, again, is de- 
scribed by 12 parameters, namely, 6 U-spin amplitudes, 5 
relative strong phases between them, and the weak phase 
7 which is the only common parameter among four pa- 
rameter sets. 

All 8 B^ — > P^P^ decays have actually been observed 
and their branching ratios and CP asymmetries have 
been measured, though, with the present statistics in 
most cases the errors are rather large. This is especially 



so for the CP-asymmetries. In any case, with 16 data 
points and 12 fit parameters one can perform a fit and 
extract the preferred values for all parameters. 

In the other 3 subsets some modes have not yet been 
observed but upper limits on their branching ratios were 
reported. Needless to say, direct CP asymmetries for 
these modes have not been determined yet. For some of 
these modes a central value and a large uncertainty are 
known. For the others, where only an upper limit at 90% 
confidence level is reported, one can take central value as 
equal to and approximately estimate the uncertainty by 
dividing the upper limit value by 2. For example, from 
B{B^ (jjp'^) < 16 we crudely estimate that B{ujp'^) = 
0.0 ± 8.0 [13. The data from upper limits helps in two 
ways. First of all, it provides additional data points, 
making a U-spin fit feasible. Second, it allows us to verify 
that the resulting fit is consistent with the current upper 
limits. 

In the case of V^P^ decays, for instance, 6 out of 8 
modes have been observed and provide 12 data points. 

The remaining two decays, K and ^tt"*", have not 

yet been observed. At present only the upper limits 
for these two modes are known: B{B+ K K+) = 



+1.3+0.6 
-0.0-0.0 



(< 5.3) 14] and B{B+ 07r+) < 0.41 [H 



0.0 

From these measurements we can crudely estimate that 



B{B+ 
0.0 ± 0.2. 



ir*°i^+) 



0.00lJ;;^g and B{B+ 



To make sure that the fit is consistent with 
the upper limits on the K and 07r+ branching ra- 

tios we add these two data points to the fit. Thus, the 
12-parameter V°P+ U -spin fit features 14 data points, 
making 7 extraction possible. 

Similarly, in the y'^l^"'" sector 5 modes have been de- 
tected and the first measurement of their CP asymme- 
tries has been attempted (though, again, with rather 
large errors) for a total of 10 data points. The other 
3 modes have not yet been observed but the upper lim- 
its were reported, allowing estimates of their branching 
ratios. The total number of V^V'^ data points rises to 
13. 

The least is known about P^V'^ decays. Not even an 
upper limit is known for K^K*'^ . Of the remaining 7 de- 
cays modes only 4 have been detected, providing 8 data 
points. For the other three an estimate of the branch- 
ing ratio can be made using current upper limits. Thus, 
there are only 11 data points and a reasonable 12 pa- 
rameter U-spin fit cannot be performed. To avoid this 
problem, one can make a joint U-spin fit to two M^M'^ 
decay subsets, e.g. a fit to both V'^P^ and P^V^ decays. 
With 7 being the only common parameter for both pa- 
rameter sets, there are 11 completely free P^V'^ U-spin 
parameters (amplitudes and strong phases) that describe 
11 P^V'^ data points. There is just enough data to make 
the joint fit work. 

TableQlshows the results of the U-spin fits to four sub- 
sets of M^M'^ decays and their combinations. The top 
part of the table shows three fits to individual subsets 
(y°P+, P°P+, yOy+) and one joint fit. As was men- 
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TABLE I: Results of the U-spin fits to various subsets of cfiarmless — ^ M'^M^ decays. Tlie bottom panel shows 7 as 
determined from direct measurements in B ^ £>(*)}{(*) decays, from indirect constraints on the apex of the unitarity triangle, 
and from SU(3) fits to charmless PP decays. 



Fit 


Subset 








Modes 


xVdof 




1. 


V°P+ 








+ <;/>7r+ K''°Tv+ p°K+ ujK+ cj>K+ 


3.97/2 


(30111)° 


2. 


pOp+ 






+ 

n n r]TV 




3.01/4 


(681?!)° 


3. 


v°v+ 






p°p+ijp^ 


■ 4>p+ K*°p+ p°K*+ uK*+ (l>K*+ 


0.05/1 


(40ljf)° 


4A. 


pOy+ 






7r°p^ rjp^ 


■n'p+ K°p+ Tv°K*+ T)K''+ r]'K*+ 


insufficient data 


4B. 


{V°P+ 


UP°v+) 








4.03/2 


(30111)° 


5. 


{V°P+ 


UP°p+) 








10.02/7 


(5411?)° 


6. 


{V°P+[jP°P+[jV°V+) 






10.47/9 


(5411?)° 



Direct measurements, BaBar 
Direct measurements. Belle ^] 
Indirect constraints, CKMFitter [20] 
Indirect constraints, UTFit [^ 
SU(3) fits to VP decays [5] 
SU(3) fits to PP decays [5] 



(67±28± 13± 11)° 
(6811* ±13 ±11)° 

(57.9 ± 7.4)° 
(66.213:8 ±0.1)° 

(59±9±2)° 



tioned before, the only way to explore P^V~^ data is to 
make a joint fit, for example, the (F°P+1JP°F+) one. 

Among the first four fits, the V^°P+ one stands out. It 
is the only fit that features a deep minimum at its pre- 
ferred value of 7, with both upper and lower uncertainties 
staying under 20°. The P"P+ fit has, on the other hand, 
a shallow minimum with very large upper uncertainty. Of 
the other two U-spin fits in the top part of the table, the 
V^V^ one produces a very shallow minimum, leaving 7 
practically undetermined. As was mentioned before, the 
current data on P°V^^ is insufficient for a U-spin fit. In- 
stead, subset was combined with the other VP 
decays (y°P+), making the joint fit possible. However, 
the effect of the data appears to be insignificant; 
the joint fit produces practically identical results to those 
of the y°P+ fit. One can draw the conclusion that the 
quality of the current data for P'~'V~^ and V^V^ is un- 
likely to significantly affect joint fits. 

This is confirmed in the lower part of the table. The 
best U-spin fit is achieved when V^P^ and P^P^ data 
are combined (30 data points) and fitted with 23 pa- 
rameters (two sets of six U-spin amplitudes and five 
strong phases, plus the weak phase 7). The joint 
{yOp+ y p^p+) fit prefers a value of 7 that is in between 
the values favored by V^P^ and P^P^ fits, namely, 
7 = (5411?)°. The addition of the V^V+ subset does 
not change this result, as expected. 

The above results are based on the latest world aver- 
age values for branching ratios and CP asymmetries in 
charged charmless B decays When the individual 

values from BaBar and Belle are very different, we em- 
ployed the PDG scaling factor S to boost uncertainties on 
the weighted averages. This modification only slightly af- 
fects the final result. The joint U-spin (V^°P+ IJ P''P+) 
fit to the unsealed data prefers the same central value 



but slightly smaller uncertainties for the weak phase: 
7=(54llJ)°. 

We also explored the joint (F°P+ IJ P"P+) fit in some 
detail with the purpose of estimating the expected im- 
provement of 7 extraction as higher statistics on B de- 
cays get accumulated. We tried to identify some specific 
modes where smaller uncertainties on branching ratios 
would help reduce the error on 7. 

We found that setting a stricter upper limit on the (f)n^ 
branching ratio is of particular importance. The current 
upper hmit B{B+ (j)iT+) < 0.41 [H is based on 89 
million BB pairs. Both B factories will each accumu- 
late in excess of 500 million BB pairs by the summer 
of 2006. The available statistics on S(07r+) decays will 
increase by about a factor of 10, leading to uncertainties 
that are about 3 times smaller than the current ones. 
With the new data point of S(07r+) = 0.00 ± 0.07 the 
joint (I/°P+ U P^P^) U-spin fit features a rather deep 
minimum with uncertainties on 7 at the level of 8°. The 

improvements in 7 extraction due to stricter upper limits 

— *o I 

on K are somewhat smaller. 

Finally, we scaled down uncertainties in all data points 
to the levels corresponding to 1 billion BB pairs. The U- 
spin fit becomes much deeper at its minimum and 7 gets 
extracted with a 6° uncertainty. Theoretical uncertain- 
ties associated with this method are expected to be small 
so it has the potential to put rather stringent constraints 
on the weak phase 7 [T7l |. 

Summarizing, with current statistics, the best U-spin 
fits allow the extraction of 7 with a reasonable accu- 
racy and the preferred value is 7 = (54ll?)° which is 
quite consistent with the current indirect determinations 
that expect 7 to lie between 42° and 73° IMEll. Note 
that the intrinsic theoretical uncertainty associated with 
possible U-spin breaking effects is expected to be rather 
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small and that U-spin symmetry is the only assumption 
that is made in this approach |22l| . Clearly, as data with 
higher statistics becomes available, the statistical uncer- 
tainties on 7 will become even smaller. At the moment 
the difference between the four values of 7 extracted from 
the four subsets is not very meaningful due to large un- 
certainties (Table nj. When all branching ratios and CP 
asymmetries in charged B decays are experimentally de- 
termined with high accuracy, U-spin approach should en- 
able extraction of 7 quite precisely from each of the four 
subsets of data. The resulting spread in 7 values should 
be small and could perhaps be used to indicate the sys- 
tematic errors inherent in the method due to residual 
U-spin breaking effects. The crucial advantage of the 



method is that the extraction of 7 is completely model 
independent and entirely data driven. Note also that un- 
like the use of isospin for a, electroweak penguins are not 
a problem in our approach. Penguin contributions are en- 
tering in an important way in this U-spin approach for 
getting 7. That means that this method is sensitive to 
new physics in the loops. In contrast, recall that the stan- 
dard B DK methods 1] involve only tree B decays. 
Comparison of 7 from these two methods is therefore im- 
portant for uncovering new physics. 

We thank J. G. Smith for helpful discussions. This 
research was supported in part by DOE contract Nos.DE- 
FG02-04ER41291(BNL). 
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